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Abstract 

We prove the Central Limit Theorem for the number of eigenval- 
ues near the spectrum edge for hermitian ensembles of random matri- 
ces. To derive our results, we use a general theorem, essentially due 
to Costin and Lebowitz, concerning the Gaussian fluctuation of the 
number of particles in random point fields with determinantal correla- 
tion functions. As another corollary of Costin-Lebowitz Theorem we 
prove CLT for the empirical distribution function of the eigenvalues 
of random matrices from classical compact groups. 

1 Introduction and Formulation of Results 

Random hermitian matrices were introduced in mathematical physics by 
Wigner in the fifties ([Wigl], [Wig2]). The main motivation of pioneers 
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in this field was to obtain a better understanding of the statistical behavior 
of energy levels of heavy nuclei. An archetypical example of random matri- 
ces is the Gaussian Unitary Ensemble (G.U.E.) which can be defined by the 
probability distribution on a space of n-dimensional hermitian matrices as 

P(dA) = const n • e~ 2n Tracc A 'dA. (1.1) 

Here dA is the Lebesgue measure on n 2 -parameters set 

{Re dij, 1 < i < j < n; Im a^, 1 < i < j < n, an, 1 < i < n} (1.2) 

2 

and const n = (7rn) - ^" • 2 n ^ 2n ~ 1 ^ 2 is a normalization constant. (1.1) im- 
plies that matrix entries (1.2) are independent Gaussian random variables 
N(0, 1 g^' 3 ) . It is well known that the G.U.E. is the only ensemble of her- 
mitian random matrices (up to a trivial rescaling) that satisfies both of the 
following properties: 

(1) probability distribution P(dA) is invariant under unitary transforma- 
tion 

A -> U^AU, U E U(n), 

(2) matrix entries up from the diagonal are independent random variables 
(see [Me], Ch. 2). 

The n eigenvalues, all real, of hermitian matrix A will be denoted by Ai, 
A 2 , . . • , A n . For the formulas for their joint distribution density p n (Xi, . . . , A n ) 
and /c-point correlation functions p n ,fc(Ai, . . . , A&) we refer to Mehta's book 
[Me]. One has 

n 

p„(Ai, . . . , A n ) = const^ • |[ |A i -A i | 2 -ea;p(-2n-^A 2 ) (1.3) 



nl I 

Pn,fc(Al, . . . , A fc ) : = — / p n (Ai, • • • , \ n )d\ k+ i . . . d\ n 

(n-k)\J R n-k (1.4) 

= det(K n (\ u \ 3 ))l J=1 , 
where K n (x,y) is a projection kernel, 

rt-l 
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and 



ibAx) = — j r • exp — • -r- exp(— x ) (1.6) 

tt3 • (2<-^!)3 V 2 / dx* 

i = 0, 1, ... , are Weber- Hermite functions. The global behavior of eigen- 
values is governed by the celebrated semicircle law, which states that the 
empirical distribution function of the eigenvalues weakly converges to a non- 
random (Wigner) distribution: 

^n(A) = -#{A, < A} JF(A) = [ X p( x )dx (1.7) 

with probability one ([Wigl], [Wig2]) where the spectral density p is given 
by 

(^VT^¥, \t\<i . . 

p(t) = < n ' - 1.8 

[0, \t\ > 1. y J 

To study the local behavior of eigenvalues near an arbitrary point in the 
spectrum x G [—1, 1], one has to consider rescaling 

\ l = x + ^l— i = l,...,k (1.9) 

Pn,l( X ) 

and study the rescaled /c-point correlation functions 

Rn.kiVl, ■ ■ ■ , Vk) •■= (Pn,l( X ))~ k ■ Pn,fc(Al, . . . , A fe ). (1.10) 

The biggest interest is paid to the asymptotics of rescaled correlation func- 
tions when n goes to infinity. For G.U.E. the answer can be obtained from 
the Plancherel-Rotach asymptotic formulas for Hermite polynomials ([PR]): 

lim R n , k (yi, ■■■ ,Vk)= Pk(yi, ■■■ ,Vk) = det(K(y i} yj)) k i>j=l . (1.11) 

The K actually also depends on x but in a very simple way. It can be 
represented as 

^^M^wdl, (L12) 
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where for all |x| < 1 the function A is just sm ^ , and for x — ±1 it is 

Ai(±y) = \ J q cos Q* 3 ± y*) d *- ( L13 ) 

The function defined by (1.13) is known as the Airy function and the kernel 
(1.12)-(1.13) is known as the Airy kernel (see [Me], [TW1], [F]). The limiting 
correlation functions (1.11)-(1.13) determine a random point field on the real 
line, i.e., probability measure on the Borel cr-algebra of the space of locally 
finite configurations, 

n = {u: = ( yi ) iez : VT > #{ Vi : \ Vi \ <T}< oo}. (1.14) 

The distribution of random point field is uniquely defined by the generating 
function 

k 

ip(z 1 , ... ,z k ; h, . . . , 4) = E Y[ z]\ 

3=1 

where Ij, j = 1, . . . ,k, are disjoint intervals on the real line, Vj = #{yi G 
Ij} = #(ij), the number of particles in Ij, and k E Z\. It follows from the 
general theory of existence and uniqueness for random point fields ([LI], [L2]), 
that if K(y, z) is locally bounded than determinantal correlation functions 
uniquelly determine random point field assuming that such randompoint fiels 
exists. The generating function ip(zi, . . . ,Zk) is given by Fredholm determi- 
nant of the intergal operator in L 2 (M 1 ): 

k 

tJ>(z!,... ,z k ) = det(5(x - y) + ^(^- - 1) • K(x,y) ■ Xijiv)), (!- 15 ) 

where xij is an indicator of Ij. 

In particular these results are applicable to the Airy kernel (1.12)— (1.13). 
We shall call the corresponding random point field the Airy random point 
field. For one-level density formulas (1. 11)— (1.13) produce 

Pi(y) = -y(Ai) 2 (y) + (Ai'(y)) 2 . (1.16) 

The asymptotic expansion of the Airy function is well known (see [01]). One 
can deduce from it 
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r • exp(— 4yz / 6) as y — > +oo. 
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Pi satisfies the third order differential equation 



p'; , (y) = -2pi(y)+4yp' 1 (y). (1.18) 

One can think about the one-point correlation function as a level density, 
since for any interval / C R 1 we have E#(I) = JjPi{y)dy. It follows from 

(1.17) that £#((-T, +00)) is finite for any T and £#((-T, +00)) ~ ^ + 

3 

0(1) when T goes to +00. The last formula means that E#((—T, +00)) — ^^- 
stays bounded for large positive T. Let us denote v\(T) := #{yi > —T} = 
#((-T, +00)), u k (T) := —(k - 1)T])), k = 2, 3, ... . Theorem 1 

establishes the Central Limit Theorem for is k (T). 

Theorem 1 The variance of u k {T) grows logarithmically 

Vaxi/ fc (T)~^2-logT + 0(l), 
and the sequence of normalized random variables Uk (p~ Euk SD- converqes in 

^/Var v k (T) 

distribution to the centalized gaussian random sequence with the covari- 
ance function E£ k & = 8 k ,i - 1/2 5 k j + i - 1/2 

Remark 1. The first result about Gaussian fluctuation of the number of 
particles in random matrix model was established by Costin and Lebowitz 
([CL]) for the kernel ""g^ . See §2 for a more detailed discussion. 
Remark 2. Basor and Widom ([BaW]) recently proves the Central Limit 
Theorem for a large class of smooth linear statistics X^J^oo f(Vi/T) where / 
satisfies some decay and differentiality conditions. Similar results for smooth 
linear statistics in other random matrix ensembles were proven in [Sp],[DS], 
[Jol], [SiSol], [KKP], [SiSo2], [Ba], [BF], [BdMK] , [Br]; see also [Sol] for the 
results about global distribution of spacings. 

Another class of random hermitian matrices, called Laguerre ensemble, 
was introduced by Bronk in [Br]. This one is the ensemble of positive n x n 
hermitian matrices. Any positive hermitian matrix H can be represented 
as H — AA*, where A is some complex valued n x n matrix and A* is its 
conjugate. The distribution on such matrices is defined as 

P{dH) = const" • exp(-n • Trace A- A*)- [det(AA*)} a dA, (1.19) 

where a > — 1 and dA is Lebesgue measure on 2n 2 -dimensional space of 
complex matrices. The joint distribution of n (positive) eigenvalues of H is 
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given by 



Pn (A 1 ,...,A n )=consC-exp -n-^A, -J]Af- J] (A, - A,) 2 . 

\ i=l / i=l l<«<j<n 

(1.20) 

The Vandermonde factor in (1.20) implies that correlation functions still have 
the determinantal form (1.4) with the kernel 

n-l 



K n (x,y) = n ■ ^2<f) e (nx) ■ </> e (ny), (1.21) 



where the sequence {<pe(x)} is obtained by orthonormalizing the sequence 
{x k ■ X2 ■ e~ x / 2 } on (0, +oo). The limiting level density is supported on [0, 1] 
and given by the equation 

Pl (x) = ^-x-*>{l-x)*. (1.22) 

(It is not surprising to see that (1.22) is the density of a square of the Wigner 
random variable!) Plancherel-Rotach type asymptotics for Laguerre polyno- 
mials ([E], [PR]) imply that by scaling the kernel K n (x,y) in the bulk of 
the spectrum, we obtain the sine kernel ^t^z^ > an d by scaling at x — 1 
("soft edge"), we obtain the Airy kernel. As we already know, the same 
kernels appear after rescaling in G.U.E. This feature, called universality of 
local correlations, has been established recently for a variety of ensembles 
of hermitian random matrices (see [PS],[BI],[DKMVZ],[Jo2],[So2],[BZ]). To 
scale the kernel at the "hard edge" x = 0, we need an asymptotic formula of 
Hilb's type (see [E]), which leads to 

lim —K (— — \ = jgfaffi ! y^L ' J '^^ ~ yftJ'M " J °(^) (i 23) 
«™ 4n n \An An) 2(x - y) ' 

where J a is the Bessel function of order a ([F], [TW2], [Ba]). The kernel 
(1.23) is also known to appear at hard edges in the Jacobi ensemble ([NW]). 
For a quick reference, we note that in the Jacobi case, a sequence {(f>e(x)} 
from (1.21) is obtained by orthonormalizing {x h (l — + x) 2 }. The 

random point field on [0, +00] with the determinantal correlation functions 
defined by (1.23) will be referred to as the Bessel random point field. There 



6 



is a general belief among people working in random matrix theory that in 
the same way as the sine kernel appears to be a universal limit in the bulk 
of the spectrum for random hermitian matrices, Airy and Bessel kernels are 
universal limits at the soft and hard edge of the spectrum. The next theorem 
establishes the CLT for v k {T) = #(((*; - 1)T, kT\), k = 1,2, 

Theorem 2 Let v{T) be the number of particles in (0, T) for the Bessel 
random point field. Then 

E v k (T) ~ -T^(k 1 ' 2 - (k - 1) 1/2 + 0(1), 

7T 

Var z/ fe (T)~^logT + 0(l), 
and the sequence of the the normalized random variable Uk ^)~ E u k ^l con _ 

^Var u k {T) 

verges in distribution to the gaussian random sequence from the Theorem 1. 

Theorems 1 and 2, as well as similar results for the random fields arising 
from the classical compact groups (see [Sol]) are the corollaries of the general 
result about determinantal random point fields, which is essentially due to 
Costin and Lebowitz. Recently a number of discrete determinantal random 
point fields appeared in two-dimensional growth models ( [Jo3],[Jo5]), asymp- 
totics of Plancherel measures on symmetric groups and the representation 
theory of the infinite symmetric group ([BOl], [B02],[BOO],[Jo4],[Okl],[Ok2]) 
If one can show the infinite growth of the variance of the number of parti- 
cles in these models ( the goal which may be probably attainable since the 
asymptotics of the discrete orthogonal polynomials arising in some of these 
problems are known) the Costin- Lebowitz theorem should work there as 
well. 

The rest of the paper is organized as follows. We discuss the general 
(Costin-Lebowitz) theorem in §2. Theorems 1 and 2 will be proven in §3 
and 54. In the Bessel case, we will see that the kernel ± sin^+y) 

naturally appears in our considerations. We recall in §4 that the sine kernel 
also appears in the limiting distribution of eigenvalues in unitary group and 
the even and odd sine kernels appear in the distribution of eigenvalues in 
orthogonal and symplectic groups and then prove Theorems 3-6 the Gaussian 
fluctuation for the number of eigenvalues in these models in Theorem 3-6. 

It is a pleasure to thank Prof. Ya. Sinai, who drew my attention to the 
preprint of Costin-Lebowitz paper some time ago, and the organizers of the 
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Special Session on Integrable Systems and Random Matrix Theory (AMS 
Meeting at Tucson, November 13-15, 1998) and the Introductory Workshop 
on Random Matrix Models (MSRI, Berkeley, January 19-23, 1999) for the 
opportunity to attend the meetings, where the idea of this paper has been 
finalized. The work was partially supported by the Euler stipend from the 
German Mathematical Society. 

2 The Central Limit Theorem for Determi- 
nantal Random Point Fields 

Let {Vt}^ 1 be a family of random point fields on the real line such that 
their correlation functions have determinantal form at the r.h.s. of (1.11) 
with kernels K t (y,z), and {It}t£R\ a se ^ °f intervals. We denote by A t an 
integral operator on I t with the kernel K t (y,z), A t : L 2 (I t ) — > L 2 (I t ), by u t 
the number of particles in I t , u t = #(/*), and by E t , Var t the mathematical 
expectation and variance with respect to the probability distribution of the 
random field Vt- In many applications the random point field Vt, and there- 
fore the kernel K t will be the same for all t. In such situations the interval 
I t will be expanding. 

Theorem (O. Costin, J. Lebowitz) Let A t = K t -Xi t be a family of trace class 
operators associated with determinantal random point fields {V t } such that 
Va.it v t = Trace(A t — A 2 ) goes to infinity as t — > +oo. Then the distribution 
of the normalized random variable with respect to the random point 

field V t weakly converges to the normal law N(0, 1). 

Remark 3. The result has been proven by Costin and Lebowitz when K t (x, y) = 
S ™(l-y)' 1 f° r an y * anc ^ 1-^*1 \ — 00 ( see IpL])- The original paper contains a 
remark, due to Widom, that the result holds for more general kernels. 
Remark 4- There is a general result that a (locally) trace class operator K 
defines a determinantal random point field iff < K < 1 (see [ So4] or, for 
a slightly weaker version, [Ma]). 

The idea of the proof is very clear and consists of two parts. Let us denote 
the £th cumulant of u t by Cg(i/ t ). We remind that by definition 



J2CSzY/£\=\og(Etexp( iZ Ut)) 
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Lemma 1 The following recursive relation holds for any £ > 2: 

e-i 

C t (y t ) = (-I)* ■ (£ - l)!Trace(A - A l t ) + £ a se C s (u t ), (2.1) 



s=2 



where a S £, 2 < s < £ — 1, are some combinatorial coefficients (irrelevant for 
our purposes). 

The proof can be found in [CL] or [Sol], §2; (of course one has to replace 
everywhere S1 °^~^ by K t (x, y)). For the convinience of the reader we sketch 
the main ideas here. We start by introducing the Ursell (cluster) functions : 

n(xi) = p(xi}, r 2 (xi,x 2 ) = P2(xi,x 2 ) - p 1 (x 1 )p(x 2 ), 

and, in general, 

k m 

^....x^^Q-rV-l)! J] r Gi {x{G s )) (2.2) 

m=l G j=l 

where G is a partition of indices {1,2,... , k} into m subgroups 
Gi,...G m , and x(Gj) stands for the collection of Xi with indices 
in Gj . 

It appears that the integral of k-point Ursell function r^Xi, . . . ,x^) 
over /c-dimensional cube I t x . . . I t is equal to the linear combination of 
Cj(v t ), j — 1, . . . k.. Namely, let us denote 

T k (v t )= ■■■ r k (x 1 ,...,x k )dx 1 ...dx k 
Jit Jit 

Then 

oo oo 

C k {iz) k /k\ = ^(exp(^) - l) k T k {v t )/k\ (2.3) 

k k=l 

Taking into account that for the determinantal random point fields 

T k (u t ) = (-l) k -(k-l)\Trace(A t ) k , 

the last two equations imply (2.1). The next lemma allows us to estimate 
Trace (A t - A\). 
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Lemma 2 < Trace(A - A\) < {£ - 1) • Trace (A t - A 2 t ). 

The proof is elementary: < Trace (A t -A e t ) = Y?~Jy Trace (A j t -A> t +1 ) < 

£ ■=} Pr'll • Trace (A t - Af) <(£-!)■ Trace (A t - A 2 t ). □ 
As a corollary of the lemmas we have Cg(i> t ) = 0(C2(ft)) for any £ > 2. 
Since Ciiyf) = Trace (A t — A 2 ) — > +oo, we conclude that for £ > 2, 

t— »oo 

fi I vt-Evt \ _ Ce(vt) n 

U H V / U ^J ((C 2 (i/ t ))V2 U - 

At the same time the first two cumulants of the normalized random vari- 
able are and 1, respectively. The convergence of cumulants implies the con- 
vergence of moments to the moments of N(0, 1). The theorem is proven. □ 

To generalize the Costin-Lebowitz theorem to the case of several intervals 
we consider l[ m \ m + 1, . . . , s, disjoint intervals of the real line, and define 
jyW = #(/M). The equation 

Y^ C ^-^i) kl /ki\ ■ ■ ■ {iz s ) k °/k s \ = log (^exp^z/W + . . . + z s v\ a) )j) 

(2.4) 

defines the joint cumulants of z/f^'s. 

Proposition 1 Ck u ...,k a i^t~\--- ^ u t(s)^ is equal to the linear combination 
of the traces 

TraceK t . X ^ ) .K t . X ^ ) ...K t . X ( j; ) 

with some combinatorial coefficients (irrelevant for our purposes), such that 
for any kj that is greater than zero at least one indicator in each term of the 
linear combination is the indicator of 

The proof immedeately follows from the analogue of (2.4) for the case of 
a several intervals. 

In the next section we will apply theses results to prove Theorem 1. 

3 Proof of Theorem 1 

For the most part of the section we will study the case of one interval 
(— T, +00). We start by recalling the asymptotic expansion of Airy func- 
tion for large positive and negative y (see [01]). 

— TTZ 00 

27T2.|y|4 s=Q z 
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\y\* 



2 ^ s=0 Z 



1 ( °° 

AHi/l) - -^-r • - (™ + 1) ■ J2(-iy ■ u f s 

7T2 • \y\4 I \ 4/ ^ Z 

OO N 



(3.3) 



sin 7T2 



III/' OO 

^J(- M) ~J2^. + 1). 

00 

+ i)-E(- 1 )' +1 '" 2 * +1 



— COS 7TZ 

I 

s=0 



y2s + l 



(3.4) 



where z = £y - \y\*; Uq = v = 1, and w s = g£± ^§S^r g3 1} , 



6s+l 
6s- 



]--u s , s > 1. In particular, as a consequence of (3.1)-(3.4) one has (1.17). 
It follows from (3.1)-(3.4) together with the boundedness of Ai(y), A[(y) on 
any compact set that for any fixed a £ I 1 all moments of #((a, +00)) are 
finite. Therefore it is enough to establish the CLT for #((— T, a)). We choose 
a = — (^ L ) 1 (y = — (^f) 5 corresponds to z = —1). We are going to show 
that the conditions of the theorem from §2 are satisfied by K ■ X(-T,a), where 
as above, this notation is reserved for the integral operator with the kernel 
K(x,y) ■ X(-T, a )(y)- 

Lemma 3 < K • X(-T,a) < 1 an d K ■ X(-t» is trace class. □ 
The kernel K(y 1: y 2 ) = AfaO-^^'fe^faO-Afe) wag obtained from 



K n (xi,X2) = v2n-^2 e=Q ^(y2nxi) -tpeW 2n-x 2 ) after rescaling Xi = 1 + -^-, 
% — 1,2, and taking the limit n — > 00. The convergence is uniform on compact 
sets. As a projection operation K n satisfies < K n < 1. We immediately 
conclude that K ■ X(-T,a) satisfies the same inequalities and since the kernel 
is continuous and non-negative definite the operator is trace class (see e.g 
[GK] or [RS], section XI. 4). Now the main step of the proof consists of 

Proposition 2. 

Var|#|-T,-(^) 1 1 ~_lLlogT + 0(l). 
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Proof. We introduce the change of variables 

2 i 

Zi = —Vi ■ \Vi\ 2 (3.5) 

and agree to use the notations Q(zi, z 2 ), qk{z\, ■ ■ ■ , Zk), k = 1,2, ... , for 
the kernel and fc-point correlation function of the new random point field 
obtained by (3.5). It follows from (1.17) that qi (z) ~ 1 + ^es. + o(^ 2 ) for 
z — > —00, so we see that the configuration (,%) is equally spaced at —00. 
The kernel Q(zi, z 2 ) is defined by 

7T 

<5(^i, 22) = — 1 r • K (vu 2/2) 

71 AM ■ ^(y 2 ) - Aiiyi) ■ Ai(y 2 ) 



1 2/x 1 4 - 1 2/2 1 4 1/1 V2 

Formulas (3.3)-(3.4) allow us to represent Q as the sum of six kernels i = 
1, . . . , 6, with the known asymptotic expansion: Q(z±, z 2 ) = Y^=i Q^H z ii Zc i)i 
where 



11 f 00 

QV>(z u Z2)~-- j • sinvr^ - z 2 ) ■ I £ (-1)™+" 

% 3 - ^ m ' n=0 (3.6) 

/ -2m- 1 _ 2n _ 2n -2m- 1 

• u 2rn ■ v 2n ■ I z 1 • z 2 + z 1 • z 2 



1 ! 

3^ 



Q (2) (-21, Z 2 ) ~ — • -J J • COS 7r(zi + 2 2 ) 



Z 3 - 2 2 3 



" 2 " 1- ^ -" ( \ . z -2n \ 1 . 



1 (-ir^ . U2m ■ v 2n ■ U 2n ■ z - 2m ~* - n 

^m,n=0 

Q (3 \zi, z 2 ) ~ ^- • -3— — 2 • 2 cos ^tt^i + ^ • cos (ttz 2 + ^ 



(3.7) 



Zi 3 - zl 

I / i\m+n+l ~ 2m ~3 _-2n-l _-2n-l r " 2m "3 

S l — X J M 2m ^2n+l I ^1 ' Z 2 — 2 1 ' Z 2 

^m,n=0 



(3.8) 
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Q ( - 4) (Z 1 ,Z 2 ) ~ ^ • J ' 2sil1 (* Z 1 + ^) ' Sil1 (^2 + ^) • 

{ E (-D 

^m,n=0 



2 i 2 2 



4 4 
1( / ~ 2m ~3 r -2n -2ra - 2m ~3 

«2ra+l V2n \ z \ ' z 2 ~~ Z l ' Z 2 



(3.9) 



Q (5) (Z 1 ,Z 2 ) ~ — • -2~ 2 ' S in7r (^l - ^2) 

^ - zl 



f 00 

{ E (-D 



4 4 

m+n+1 I ~ 2m ~3 _-2n-l , _-2ra-l ~ 2m ~3 

• «2m+l • ^2n+l " I ^1 " ^2 + ^1 ' Z 2 



(3.10) 



Q {6) ( zi,^) ~ ^ • —2- — 2 ' C0S7r(zi + z 2 ) 
z? — z$ 



( / 4 4 \ ~» 

^m,n=0 ^ ' J 

(3.11) 



We denote by Qm, n (zi, z 2 ) the (m,n) th term in the asymptotic expansion of 
Q {i) (z u z 2 ). Then 



(1) _ sinTr^x-^) 1 : 

Vo,OV^l> ^2j - 2 2 3^ 1 1 



1 

3 

2 



1 1 



sin 7r(zi — z 2 ) zf + z-p ■ z% + z 2 3 

7r(^l -22) oil 

We note that near the diagonal Q^o i s essentially a sine kernel. We also will 
need 



2 1 1 

Qq,q{ z 1i Z< 2) — 



_ cos7r(zi + z 2 ) zl - zl ■ zl + z\ 



1 1 

3 • zl ■ z^ 
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Let us define S(z 1 ,z 2 ) = Qo,b( z u z 2 ) + Q ,b( z ii ^ 2 ), U(z 1 ,z 2 ) = Q(z 1 ,z 2 ) - 
S{z u z 2 ). 

Lemma 4. 

I l I-l (0^))V^ 2 = L-^logL + 0(l). (3.12) 



Proof. The integral can be written as 



I f f L f siring -z 2 ) 



9 Ji Ji V n(zi - z 2 ) 



2 11 2\ 2 

+ z 3 ■ ^ 2 3 + 2 2 3 
1 1 
zl ■ zl 



dz\dz 2 



2 /■ L_1 /sinvra^ 2 ' 



7TM 



9 j 

We represent the inner integral as 



Z + U \ 3 „ / Z \ 3 
' + 1 + . 

z -\- u 



dzdu. 



r fi±iiy +2 .fi±!i)* + f_!_) l +2 . 

A V * / \ * / \* + «/ v^ + w/ 



+ 3 



7i(«) + / 2 (u) + 3-(L-u-l), 



(3.13) 



with 



L—u 



h(u) = 

J 

/,(«) = ( 



2 



2 + M 



+ 2 



+ 2 



z 



z 



z + u 



dz, 



dz. 



To calculate h(u) we introduce the change of variables t = (£±^)3. Then 
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1 

(1+u) 3 



U 



= J L i(t 2 + 2t)- [y^J dt= {(l + u)3 + 2- (l + u)3 



("I) 



L \ 3 



L — u 



2 • 



L \ 3 



L — u 



-L + u) 



/•(!+«) 3 j 



(3.14) 



We have 



M 



(i+«)3 1 

1 (2t + 2) • rrft = u • 



4m , 
= — - log 

3 & 

2m, 
log 

3 & 

2m , 
+ ylog 



t 3 - 1 



(1 + M)3 - 1 



(l+u)3 4 
L I 



1 t+J 

3 U-l t 2 + t + l 



4m , 
— log 

3 & 



L \ 3 



L — u 



- 1 



(1 +M)3 + (1 + u )3 + 1 



L \ 3 



L — u 



+ 



L \ 3 



L — u 



+ 1 



and 



(3.15) 



^(u) = _(i + u ) 3 -2(1 + m)3 +L- 1 - - +2L- 1 - 



i 

M\3 



2 

«\ 3 



4m , 
+ T los 



2m , 



(1+M)3 - 1 







4m , 




-ylog 









L \ 3 



L — u 



+ 



L \ 3 



L — u 



+ 1 



(3.16) 



2m , 
ylog 



2 1 
(1 + M)3 + (1 + U )3 + 1 



In a similar way in order to calculate 
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h(u) 



1+u 



2 1 

+ 2 • ( ) dz, 



z 



we consider the change of variables t 
1 

'(i=u) 



z 



;^)s. Then 



+ 



i 



2 • I — V) • L — ((l + u)-l+2- (l + u)"i )•(! + «) 



l 

L \ ~3 



(3.17) 



/ V / (2t + 2) • -dt = -(1 + u) 3 -2(l + u) 

i(i+«) 3 r — 1 



i 



+£ . ( i_«)S + 2L .(i_«)* + | u . l0| 



L - w\ 3 



-u • log 

3 & 



1 - (1+U)"3 



2 . 
-u • log 

3 & 



L - u\ 3 



'L-1i\3 

- 1 — ^- + 1 



+ g« • log 



2 1 
(1 + U) 3 + (1 + U ) 3 + 1 



It follows from (3.13) that 



L J-L 



Qofi(zi,z 2 )^ dzidz 2 



L-l 



9 Jo 
We note that 



/sm7ra\ ./ 7i / u ) + /2 / u ) + 3 L _ 3u _3) du- 
\ nu J 

L-l / • \ 2 

/ sin 7TU \ 



(3.18) 



I 



, Ldu = - + 0(1); 
\ nu J 2 



(3.19) 



pL-l 
JO 



(sin iru)' 



du 



nu 



2tt 2 



logL + 0(l); 



(3.20) 
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lu ( u 



du 



L 
2 



(3.21) 



6tt 2 



logL + 0(l); 



L-l 







/ sin 7tu \ 
\ iru ) 



■ L ■ 



(l-£) 5 d«=£--^log£ + 0(l). (3.22) 



2 3tt 2 



The combined contribution of all other terms to (3.18) is 0(1). Indeed, 

" ; 1 f sin n \ 2 



•(1+^)3 =0(1); 



f L ~ x /sinTraY 
Jo \ nu J 
f 1 ' 1 /sinvnA 2 . i 



Jo \ nu J 3 



(1+U)3 - 1 



2w , 
ylog 



(1 +U)3 + (1 +u)3 + 1 



/ Sin 7TM 

\ 7ra 



■)' 



2 Z^ -1 sin 2 7TU 



37T 2 J 1i 2 U 



2u f L 1 sin 2 7tu 



log 



(1+U)3 - 1 



2 1 
(1 +U)3 + (1 +u)3 + 1 



3vr 2 7 



(u"3 ) du = 0(l); 



L-l 



sin 7T-u \ 4-u 
iru J 3 



log 



1-(1 + U)"3 



du = 0(1); 



sin 7TU \ 2u . 
--5- log 

7TU / 3 



_2 1 
(1 +U) 3 + (1 + W) 3 +l 



du = 0(1). 
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The last expression to consider is 



* T 1 / \ 2 

' 'sinvruA 4 



TTU 



U ■ log 



2 , 
-u ■ log 

3 & 



L \ 3 



L — u 



+ 



L — u 
L 



- 1 



L — u 



+ 1 



du + 



L-l / • \ 2 

' Sin 7TM \ 





L-l 



TTU 



log 



1 - 



L - u\ 3 



L-l 



/ Sill 7TM \ 2 



du- I —u log 

\ TTU J 3 



sin 7ru \ 4 

• -w 

nil J 3 



L - m\ 3 



L — m\ 3 



2 

3^. /0 



L 1 sin 2 7rit 



it 



•log 



:^) 5 -i 



X 



(i^) 5 +(^) 5 +l 



3tt 2 



; 1 sin 2 7tm 



log 



L — u 



du 



fc=i 

k 1 °° 1 

fe=l 

Combining all the above integrals and looking specifically for the contribu- 
tions from (3.19)-(3.22) we obtain 

/— 1 r— 1 , > 2 1 { L 1 Z/ 

y (Qo3(*i> ^i^2 = g( 3 'Y _3 '272 logL+ 2 - 



6tt 2 



log L + 2 • - - 2 • At • log L + l — log L + 2 • - 

5 2 3tt 2 & 2 3vr 2 & 2 



-2 



6tt 2 



log L + 0(1) 



L — 



3tt 2 



log L + 0(1). 



□ 
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In the next lemma we evaluate integrals involving Q^ . 
Lemma 5 

a ) I-l 5-1 (Qo?o(*i> * 2 )) ^1^2 = ^ logL + 0(1). 

b ) 5-1 5-1 Qvl( z ^ z z) ■ Qofi( z i, z 2 )dz l dz 2 = 0(1). 
Proof. The integral in part a) can be written as 



1 

9 h 



L r L 




COS7r(^i + Z 2 ) 

ir(zi + z 2 ) 



2 1 1 2 \ 2 

ry 3 >y 3 *y 3 I ^/ 3 

Z l Z l ' *2 ~T Z 2 

1 1 

zf ■ z$ 



dz\dz 2 



1 f IL /cos7ra\ 2 r 
= 97 2 V nil ) J, 

We denote the inner integral by 



u — z 



- 1 + 



U — Z \ 3 



dzdu. 



h{u) : 



2 

2 \ 3 



U — Z 



u — z 



+ 



u — Z \ 3 



If — Z \ 3 



+ 1 



= 2 • 



U — Z\ 3 



- 2 



It — Z \ 3 



dz + (u- 1) 



2 'Li (i2 - 2t) '(^)' <ii+(i, - i) 
i 1 



G "- / (/ 3 . 1), ■ { " h 



= «(i + e/ 



00 t 4 - 2t 3 



(ft +0(«3). 



(t 3 + l) 2 ' 

43 2 

Taking into account J °° jt^ydt = 0, we obtain ^(m) = w + 0(w3), and 



2L 



COS7TU, 2 2 



nil 



(u + 0(u3))du 



187T 2 



logL + 0(l). 
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Now let us consider the integral in part b): 
" L rL " ~(zi + z 2 ) sin7r(2;i - z 2 ) 



i r r cosTT 

9 A 7i 



Z 1 + Z 2 ) 7T • (^1 - Z 2 ) 



1 

,3 
1 

-2L 



z 2 



Z-l 



Z 9. + 2-2 



1 

z l 



dz\dz 2 



7 3 
z 2 



i r 2L cos™ r- 2 

9 J 2 7TU J 



Sin 7Tf 



71V 



U + V 



u — V 



+ 



u 



V 



U + V 



+ 1 



dvdu. 



It is not difficult to see that oscillations of trigonometric functions make this 
integral to be of order of constant. To show this we write the inner integral 
as h(u) + 15(11) + Iq(u), where 



h{u) = 

h{u) = 
h{u) = 



u-2 



sin irv ( u + v 



71V 



U — V 



dv, 



u-2 



sin irv 



u-2 



irv 



sin irv 



irv 



U + V \ 3 



u 



dv. 



dv. 



Integration by parts gives h{u) = \ + 0(u r ). Consider now 
1 

pu— 2«4 

h(u) = / 

^0 



sin irv 



irv 



u 



u 



u-2 



dv + 



sin irv 



irv 



u 



u 



dv. 



1 

u-2«4 



(3.23) 

1 2 _j_ 

The second integral in (3.23) is 0(w4 • u~ l ■ u3) = 0(-u 12). As for the first 

one we introduce x(M even), an indicator of the set where the integer part 

of v is even, and write it in the following form 



I 



1 

u-2«4 



sin irv 



ir 



■X([v] even) 



u + v \ 3 



u 



v + 1 



u + v + 1 



u 



dv + 0(u 2 
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The absolute value of the last expression is estimated from above by 

/ 2 2 

sin7rt>| / v ■ (u — t>)3 — (u — v — 1)3 • (v + 1) 



I 



1 

«-2u4 



2 2 
U • + 1) • (ti - V) 3 • (u - V - 1) 3 



JO 



TT 

2 ( - - \ 

"""" I Sin I / (u - v)3 - [u - v - 1)3 \ 



TT 



2 2 
(U + 1) • (U - u)3 • (u - U - 1)3 



+ 



f u ~ 2u4 \f&nnv\ / 1 \ , , 1 

/ J L 2- ^ = 0(^6). 

^ 71 \V (V + 1) • (U-V)3 I 



_1 _ J_ 

Therefore i^u) = 0(it 6), and similarly I^u) = 0(u 12). As a result 



si 



sm ttu , , . r , . _ / \ \ 7 1 f sin™ , 

(J 4 (w) + M«) + h{u))du = — du 

18 J 2 TTU 



TTU 
r2L 



1 f sin7rw , J_ N , 



Lemma 5 is proven. □ 
As a result of the last two lemmas we have 



J J S\z 1 ,z 2 )dz 1 dz 2 = L-^logL + -^\ogL + 0(l) 



11 logL + 0(l). 



187T 2 

Remember that S was defined as S(zi, z 2 ) = Q^(zi,z 2 ) + Q^ (zi, z 2 ). To 
finish the proof of the CLT for #(— T, +00) we just need to show that the re- 
mainder term U (zi, z 2 ) = Q(zi, z 2 ) — S(zi, z 2 ) is insignificant in the following 
sense: 

Lemma 6 

a) J~l J~l U 2 (z l7 z 2 )dz x dz 2 = 0(1) 

b ) S-l I—L U 0*1 > z i) ' S ( z ^ z 2 )dzxdz 2 = 0(1) 

c) j: l L U(z,z)dz = 0(l) 

Proof. We shall establish part a). Parts b) and c) can be treated in a similar 
manner. Repeating the calculations of the last two lemmas, it is easy to see 
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that for any fixed indices (i, m, n) such that (i — — 2) + m + n > 0, we 
have 



/ / {Q%{z l)Z2 )) 2 d Zl dz 2 = 0(1). 

J-L J-L 

Let us now choose N to be sufficiently large and write U (zi, z 2 ) = Un(zi, z 2 ) + 
V N (z 1 , z 2 ), where 

2 6 

u N ( Zl ,z 2 ) = j2 E q$»(*i>*o+£ E Qgn(*l>*0- 

i=l 0<m+n<iV i=3 0<m+n<iV 

We see that J '_*(Un(zi, z 2 )) 2 dz\dz 2 = 0(1). Asymptotic formulas (3.6)- 
(3.11) imply that 



W^z^K^^-iz^ + z-™) 



Z\ Zn 



(3.24) 



< • - • (zf + zl) ■ (z~ 2N + z- 2N ) . 

\z\ — z 2 \ 2 \ / 

It follows from (3.24) that if we choose N > 2 then 

(V n (z 1 ,z 2 )) 2 dz 1 dz 2 = 0(1). (3.25) 



\zi-z 2 \>\ 



(The integration in (3.25) is over the subset of [1, L] x [1, L]) . Indeed, to 
estimate the integral over z 2 we write 



2 r L Z A 2 j 



\zi-z 2 \>^ \( z l - Z*) z l J ' Jl (Z 1 -Z 2 ) 2 + 1 Z 2 

2 

Integrating over Z\ we arrive at (3.25). To integrate V$ near the diagonal 
we observe that kernels Q(zi,z 2 ), Qm, n (zi,z 2 ) are bounded in [1, +oo) x 
[l,+oo); therefore there exists some const'^ such that |Viv(zi,Z2)| — const'^ 
and ff, ,i (Vn(zi, z 2 )) 2 dzidz 2 < r + °° 2 c °n st jv ^ z _ o(l). Lemma 6 is 

JJ \Z\— Z2|< _ 2~ 1 2 2 — 

2 2 

proven. □ 
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3 

Taking L = -^Tz we deduce from Lemmas 4-6 and (3.5) that 

2 

2 _ / 37r \ 3 

\- I Hu>e | / i '^-) ')) = f 2 K(y,y)dy 

-T 

2 2 
'3£\3 /3vr\3 



r (f)* /-(¥)- _ w w 

K {yi,y2)dyidy 2 

J-T J-T 



2 2 22 

— Ts + 1 Ts 

3tt - v ; 3tt 

11 / 2 3\ 
+ 18^ l0g U T V +2(1) 
" logT + 0(l). 



127T 2 



This finishes the proof of Proposition 2 as well as the proof of the CLT for 

#(-T,+oo). 

In a very similar way one proves the CLT for arbitrary v k {T) = # ((—kT, —(k — 1)T)) , 
fc > 1. To prove the result for the joint distribution of {uk(T) } we note that 
the decay of K(x, y) off the diagonal implies 

Ci,iMT),^(T)) = CovK(T),^(T))= / / K 2 ( Xl ,x 2 )d Xl ,dx 2 = 

J-kT J-YT 

-Tracex[-JT,-(i-i)T) • # • X[-fcr,-(fc-i)T) • # = £2(1) if |k - 1| > 1. 
This together with 

Vax^i/,(T)j =ll/(127r 2 )logT+0(l), Var(z/ fe (T)) = 11/(12tt 2 ) logT+0(l), , fc = 
implies that 

Ci,iK(r), ^(T)) = Cov K(T), ^(T)) = 

ll/(127r 2 )logT + 0(l) 
for | k — 1 1 = 1 . Therefore as T — > oo . 

v k {T)-Ev k {T) ^(T)-^(T) 

^ ; A -" d fc ,i - 1/2 dfc,j-i - 1/2 d fci , + i 

VVar u k (T) y/Var 
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To take care of the joint cumulants of higher order it is enought to prove 
Lemma 7 Let at least two indices in (k±, . . . ,JC S ) are non-zero. Then 

C kl ,..., ks {vi{T), ... ,v s (T)) = O(logT) 



Proof According to Proposition 1 Lemma 7 follows from 
Lemma 8 

Tracer i^ 2 # . . . K X i s K Xh = O(logT) 

where X i s are the indicators of the intervals (—IjT, — (lj — 1)T], lj G Z\ and 
at least two intervals are disjoint . 

Proof This has been already established for s = 2. Let s > 2 . Since not 
all indices coincide by cyclicity of the trace we may assume l\ ^ l 2 ■ Now if 
h = h we can use the positivity of Xh^Xh^Xh to write 

\Tr&ce(xh K Xh K Xh K Xu ■ ■ ■ K Xh K Xh )\ < Trace (xh K Xh K Xh) ■ \\ K Xu . . . K Xls K Xh 

< Trace ( Xh K Xl2 K Xh ) 

where we used || K ||< 1, j| xij ll< 1- Since Trace (Xh^Xh-^Xh) = 
O(logT) Lemma 8 is proven when h — h ^ l<i- 
If l\ 7^ l 3 one more trick is needed. Let us denote 

Di = Xh K Xh, D 2 = K Xh K ■■■Xh K Xh- 

Then 

Trace (xh^Xh K Xh ■ ■ ■ Xh K Xh) = TracepiAi) 
< (Trace( J B 1J B*)) 1/2 (Tmce(D 2 D;)) 1/2 
(see [RS], volume I, section VI. 6). As before 

Tracep^) = Trace ( Xll Kxi 2 Kx h ) = O(logT) 

To obtain a similar bound for Trace(-Di-Di) we define 1 < p < s as the 
maximal index such that l p ^ h ■ Since we assume in Lemma 8 that there 
are at least two different indices, such p always exists. Then 

Trace (K X i 3 K . . . Xi p KXh ■ ■ ■ K Xh ■ ■ ■ K Xh)' {Kxi 3 K . . . Xi^Xh ■ ■ ■ K Xh ■ ■ ■ K Xh)* = 
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Trace (K Xh K . . . X i p K Xh ■ ■ ■ K Xh ■ ■ ■ K Xh )- {xhK . . . XhK ■ ■ ■ K Xh K Xlp . . . K Xh K) . 
Using the identity Trace(-Di_D 2 ) = Trace(_D 2 -Di) where 

Di = K Xh K . . . X i p _,K, D 2 = Xlp K Xh K . . . Xh K Xh . . . ^K^K^^K . . . K Xh K, 
we can rewrite and estimate the r.h.s. as 

(Trace ( X i p K Xh ... Xh K... Xh K Xlp ) ■ . . . K Xh KK Xh K . . . X i p _,K) \ 

< Trace { Xlp K Xh . . . K Xh K . . . Xh K Xlp ) ■ || K Xlp _ x . . . K Xh KK Xh K . . . Xlp _ x K \\ 

Here we used the positivity of X i p K X i x K X i x . . . K X h ■ ■ ■ K\h ■ The norm 
of the last factor is again not greater than 1. Finally, 

Trace { Xlp K Xh K Xh . . . K Xh . . . Xh K Xlp ) = Trace { Xh K Xlp K Xh . . . K Xh . . . K Xh ) < 

Trace ( Xh K Xlp K Xh ) ■ 1 = O(logT). 

Here we also used cyclicity of the trace. Combining the estimates for TraceDiD^ 
and Trace-Di-DJ we finish the proof of the Lemmas 7 and 8. □ 
It follows from Lemma 7 that the higher joint cumulants of the normalized 
random variables go to zero which implies that the limiting distribution func- 
tion is gaussian with the known covariance function. Theorem 1 is proven 

□ 

4 Proof of Theorem 2 and Similar 

Results for the Classical Compact Groups 

The Bessel kernel has the form (see §1) 

K(M2) = ' (4.1) 

2/i, 2/2 e (0, +oo), a > -1. 
The level density is given by 

Pl (y) = K(y,y) = l -J a {^yf - \j a+ i(Vv) • Ja-i(Vy)- (4.2) 
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The asymptotic formula for large y is well known in the case of Bessel func- 
tions (see asymptotic expansion in (4.5) below). In particular, one can see 
that 



pi(y) 



1 



for y — > +00. 



(4.3) 



The last formula suggests to make the (unfolding) change of variables Zi = 
% = 1,2. The kernel Q corresponding to the new evenly spaced random 
point field is given by 



1 1 

Q(z 1 , z 2 ) = 2ny? ■ y£ ■ K(y 1 , y 2 ) 



\ \ Jqjnzi) ■ irz 2 ■ J' a {nz 2 ) - tiz 1 ■ J' a {nzi) ■ J a {nz 2 ) 
z 1 z 2 ■ 



(4.4) 



z\ - z\ 



The asymptotic expansion of the Bessel function at infinity is due to Hankel 
(see, for example, [01]). 



( 2 \ 2 

■/«(*)"( — 
\7TZ 



COS Z 



A 2s (a) 



— sin \ z 



1 1 \ 
-aix — -71 • > ( — 1) 

2 4 i ^ v ; z 

' s=0 

- -an - -7r • > (-1) , 



s=0 



(4.5) 



where 

A (a) = 1, A 8 (a) 
Similarly, 



(4a 2 - l 2 ) ■ (4a 2 - 3 2 ) . . . (4a 2 - (2s - l) 2 ) 
s\-8 s 



(4.6) 



cos 



-sin ( z - l -arr- l -^.y{-iy.^ 
( z - a7r--7r) • jT(-l) 

\ _ / <,_n 



^2s+i(«) 



(4.7) 



The coefficients B s (a) can be obtained from (4.5)-(4.6); for example, S (o;) = 
1. 
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It follows from (4.5)-(4.7) that for a = ±| 

sinTT^-^a) sinTr^ + ^-g-l) 

Wi, ^2) = - ( s— H - ( ; 7 , (4.8) 

n(zi - z 2 ) n(z 1 + z 2 ) 

which can be further simplified as 

siri7r(,2i — z 2 ) sin7r(zi + z 2 ) 
n(zi - z 2 ) T rr(z 1 + z 2 ) 

For general values of a a small remainder term appears at the r.h.s. of (4.8). 
To write the asymptotic expansion, we represent Q(zi,z 2 ) as the sum of six 

6 

kernels: Q(zi,z 2 ) = Q^{zi, z 2 ), where 
i=i 

nzi-z 2 ) n tl Q (4.9) 



n 



{z\- 2n ■ z 2 2m + z~ x 2m ■ z\- 2n ) , 



^ 1 2 ' n,m=0 

• 5 2m («) • {z\~ 2n ■ z 2 2m - zf m ■ z\- 2n ) , 

g i3j (^i,^) 



(4.10) 



^ _ 2cos (tt* - f - f) • cos (tt, 2 - f - f ) 



00 

• ("l) n+m+1 • ^2n(«) • B 2m+1 (a) ■ [zf n ■ z 2 2m + z~ x 2m ■ z~ 2 2n ) , 

n,m=0 

(4.11) 

(1) ^ 2 i ).sl.(^- ¥ -i) 

00 

• £ (-1)"+™ • A 2n+1 (a) • B 2m (a) ■ {z^~ 2n ■ z\~ 2m + z\~ 2m ■ z^~ 2n ) , 

n,m=0 

(4.12) 
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*Vl ^2 J n>m=0 (4.13) 

„ im=0 (4.14) 

r> /_,\ / -l-2n ,-2m -2m -l-2n\ 

• -D2m+l(aJ • ^i -^2 ~ - ^2 J- 

The analysis of (4.9)-(4.14) is very similar to §3. One can see that the only 
contribution to the leading term of the variance comes from Q^ Q {zi, z 2 ) + 

Q^ (zi, z 2 ) which is exactly the r.h.s. of (4.8). It can be shown by a straight- 
forward calculation that 

L 'q$(z, z) + Q$(z, zj) dz~L + 0(1), 

Q$(zi, Z2) + Q$(*i, z 2 ))fdz 1 dz 2 ~ L - ^ log L + 0(1). 
i 

Taking into account that L = we finish the proof. □ 
The kernels sm ^ x ~y~> ± sm j(*+y) are we n k nown j n Random Matrix Theory. 

n(x— y) -K(x+y) J 

For one, they are the kernels of restrictions of the sine-kernel integral operator 
to the subspaces of even and odd functions and play an important role in 
spacings distribution in G.O.E. and G.S.E. ([Me]). They also appear as the 
kernels of limiting correlation functions in orthogonal and symplectic groups 
near A = 1 ([Sol]). Let us start with the even case. Consider the normalized 
Haar measure on SO(2n). The eigenvalues of matrix M can be arranged in 
pairs: 

exp(i6>i), exp(-^i), . . . , exp(i9 n ), exp(-i9 n ), < 9 1 ,9 2 , . . . , 9 n < ir. (4.15) 

In the rescaled coordinates near the origin Xi = (2n — 1) • % — 1, . . . , n, 
the /c-point correlation functions are equal to 

, . _ / sin7r(a;j — xj) 

R n , k[Xl , ...,x k ) = te\ {2n _ 1) . sin{7r . {Xt _ Xj)/{2n _ 1)) 

sin7r(x)i + x^ 



(2n - 1) • sin(vr(xi + x j )/(2n - 1)) 



i,j=l,...k 
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In the limit n — > oo the kernel in (4.16) becomes sm j( Xt x ^ _|_ snnr(x,+x 3 ) ^ ^ 
we consider rescaling near arbitrary < 9 < it the limiting kernel will be 
just the sine kernel. 

Let us now consider the SO{2n + 1) case. The first 2n eigenvalues of 
M G SO(2n + 1) can be arranged in pairs as in (4.15). The last one equals 
I. In the rescaled coordinates near 9 = 

_ nOj . _ 

Xi , X 1, . . . , Ti, 

7T 

the /c-point correlation functions are given by the formula 

sin7r(a;j — Xj) 



Rn,k(xi, ...,x k ) = det 



2n ■ sm(n ■ (x, — Xj)/2n) 
sin7r(xi + Xj) \ 
2n ■ sin(7r • ( Xi + Xj)/2n)) i j=1 k 



(4.17) 



In the limit n — > oo the kernel in (4.17) becomes smn ( x * x ^ _ sm A x i+ x j) _ jf 
we again consider rescaling near < < 7r, the limiting kernel appears to be 
the sine kernel. The case of symplectic group is very similar. Let M € Sp(n). 
Rescaled /c-point correlation functions are given by 



Rn,k(xi, ...,x k ) = det 



Sin 7T{Xi — Xj / 



(2n + 1) • sm(ir(xi - Xj) /{2n + 1) 
sin7r(xj + Xj) 



(4.18) 



(2n + 1) • sin(7r(^ + x j )/(2n + . ; , 

One can then deduce the following result from the Costin-Lebowitz The- 
orem. 

Theorem 3 Consider the normalized Haar measure on SO(n) or Spin). 
Let 9 G [0, 7r), 5 n be such that < 5 n < n — 9 — e for some e > and 
n-5 n — > +oo. Denote by u n the number of eigenvalues in [9, 9 + 5 n ]. We have 
Eis n = ^-5 n + 0{l), 



Var i/„ 



£log(n.<J n )+0(l) t/0>O, 
2 -^log(n-5 n ) + 0(l) z/0 = O 



and i/ie normalized random variable ^= = converges in distribution to the 
normal law N(0, 1). 
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Proof. Let K n (x,y) be the kernel in (4.16), (4.17), or (4.18). We observe 
that < K n ■ < Id as a composition of projection, Fourier transform, an- 
other projection and inverse Fourier transform. To check the asymptotics of 
Var v n is an excercise which is left to the reader. □ 

The case of several intervals is treated in a similar fasion. 
Theorem 4 Let 5 n > be such that 5 n — > 0, n8 

# ((k — l)5 n , kS n ]) . Then a sequence of normalized random variables Vk f Euk - n 

converges in distribution to the centalized gaussian sequence {£k}kLi} with the 
covariance function 



S k>l - 1/2 S k>l+1 - 1/2 6 k>l - u if k > 0, I > 
S ,i-1/V2 6 ltl , if 9 = 



Finally we discuss the unitary group U(n). 
The eigenvalues of matrix M can be written as: 

exp(i0i),exp(i0 2 ),...,exp(i0 n ).cZ,O < 6 1 ,9 2 ,...,6 n < 2n. (4.19) 

In the rescaled coordinates Xi — n ■ % — 1, . . . , n, the fc-point correlation 
functions are equal to 

„ , / sin n(xi — x^ \ 
R n:k (x u ...,x k )=det( . y ^rrr • (4.20) 



n ■ sm{7r ■ ^ - Xj)/n) y iJ=l k 



In the limit n — > oo the kernel in (4.20) becomes the sine kernel. We finish 
with the analoques of the last two theorems for U(n). 

Theorem 5 Consider the normalized Haar measure on U(n). Let 9 e 
[0, 2n) , S n be such that < 5 n < 2n — 9 — e for some e > and n- 5 n — > +oo . 
Denote by u n the number of eigenvalues in [9, 9 + S n ] . We have Ev n = ^-S n , 

Var v n = \ log(n • S n ) + 0(1) 

and the normalized random variable ^===^ converges in distribution to the 
normal law N(0, 1) . 

Theorem 6 Let 5 n > be such that 5 n — > 0, n5 n — > oo and v k<n = 
# ((k — l)8 n , kS n ]) . Then a sequence of normalized random variables Uk, p Euk - n 

yVar u k>Tl 
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converges in distribution to the centalized gaussian sequence {6fc}fcLi} with the 
covariance function 

E£k£,i = h,i - I/ 2 4,m - 1/2 5 k> i-i 

. Remark 5. Results similar to Theorems 5,6 in the regime 5 n = 5 > have 
been also established by K.Wieand ([W]). 

Remark 6. For the results about smooth linear statistics in the Classical 
Compact Groups we refer the reader to [DS], [Jol], [So3]. 
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